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Our most significant advance in 1984 was the computation of 
improved versions of the multi-component VAL solar models that include 
opacities for 17 million atomic and molecular lines computed by Kurucz 
the year before. These models now produce a consistant match to 
observations in the ultraviolet and far infrared and to the Ca H and K 
lineB. A short paper was presented at the AAS meeting in Baltimore 
by Avrett, Kurucz, and Looser, "New models of the solar temperature 
minimum and low chromosphere" (Bull. Amer. Astron. Soc., 16, 450, 1984), 
and Avrett gave a review talk at the "Chromospheric Diagnostics and 
Modeling Workshop" at Sac Peak in August. A copy of that paper is 
attached as the description of the new models for this report. We are 
preparing a large paper for the Astrophysical Journal. 

This year we have continued to work on models and spectra of 
sunspots in collaboration with Per Maltby and Olav Kjeldseth Moe of 
the Institute for Theoretical Astrophysics in Oslo. Kjeldseth Moe is 
analyzing the sunspot observations in the ultraviolet region 140-168 nm 
obtained by the NRL High Resolution Telescope and Spectrograph, while 
Maltby has made extensive photometric observations of sunspot umbrae and 
penumbrae in 10 channels covering the -wavelength region 387-3800 nm. 
Maltby visited Cambridge in August to work on the models. We are just 
completing a paper, "A new sunspot umbral model and its variation with 
the solar cycle" by Maltby, Avrett, Carlsson, Kjeldseth-Moe, Kurucz, 
and Loeser. We are continuing to work on detailed calculation of the 
spectrum and comparison with observations. The Oslo group is also 
reducing HRTS data for other solar features and they will make thoBe 
data available to us, 

Kurucz spent July through December at Los Alamos and was allowed 
to use substantial amounts of time on their Crays, CDC 7600s, and 
VAXes in a collaborative arrangement for computing atomic and molecular 
opacities for cool stars and the earth's atmosphere and for doing 
atmospheric modelling. The sun is the first tea* case, both to check 
the opacity calculations against the observed solar spectrum, which for 
ground-based spectra includes the terrestrial line absorption, and to 
check the purely theoretical model calculation against the observed 
solar energy distribution. We hope these arrangements will continue 
through 1985, but we cannot plan into the future because Los Alamos 
computers are becoming more saturated. Kurucz has not been been able 
to get access to large computers at any other laboratory. It would be 
very helpful if there were some way that NASA could provide us with 
computer time, perhaps at Ames. 

Last year we discovered that ‘the 1972 Kurucz-Feytremann atomic 
line data were not able to account for the solar ultraviolet opacity 
because or missing excited configurations of iron group elements. 
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Generally, laboratory analyses of atomic spectra ore not sufficiently 
complete to predict all the lines that can be significant in the solar 
spectrum or in computing opacity. In fact, ouch analyses are often 
quite inadequate. In order to predict configurations that hove not been 
observed, or that have been only partially observed, and to estimate 
the configuration interaction among those configurations and the known 
configurations, it is necessary to make a theoretical calculation. 

During Kurucz's visits to Los Alamos he has used Robert Cowan's Hortree- 
Fock program to build up a library of predicted Slater and transition 
integrals for 50 low configurations of the first 10 ionization stages 
of the elements up through zinc, and of the first 5 ionization stages 
for all the heavier elements. During 1984 Kurucz used these data in 
making new calculations for the atomic lines, but the work haB been 
very difficult and has gone much more slowly than anticipated; Fe II 
is finished, Fe I is almost complete, Cr I and II and Ni I and II are 
still in the least squares fitting phase. Other less abundant elements 
have not yet been started but they should play only a secondary role in 
the opacity, in any case, and will be easier to compute because most of 
them have simpler spectra. We con give some idea of the complexity by 
citing the statistics for Fe II: 



even 

odd 

number of configurations 

22 

16 

number of levels 

5723 

5198 

largest Hamiltonian matrix 

1102 

1049 

number of least squares parameters 
(many fixed at Hartree-Fock) 

729 

541 

total number of lines saved 

1,112 

,322 


Each least-squares iteration took a significant amount of time on a Cray 
and many iterations were required. In these calculations radiative, 

Stark, and van der Waals damping constants and Lande g values are 
automatically produced for each line. As soon as Fe I is completed we 
plan to compute sample spectra using both the old and new Fe line listB 
to show the improvement in the ultraviolet wavelength regions that are 
most discrepant with observation. 

We made very little progress on molecules in 1984 because of the time 
needed to complete the atomic calculations. We obtained Brault's FTS 
CN spectra but have not yet reduced them. They will provide a tremendous 
advance in our knowledge of CN and should account for many weak features 
in the solar spectrum. We got the Air Force atmospheric line list working. 
There are many weak atmospheric lines that are clearly present in our high 
quality solar spectra but that are not in the line list. We are also 
trying to get all the old Los Alamos molecular line data into shape for 
use in our opacity calculations. 

We also did the following work in 1984. 

Peter UlmBchneider asked Kurucz to produce CO infrared opacities 
for use in acoustic wave chromospheric heating calculations. We have 
computed the absorption coefficient spectra per CO molecule for a range 
of temperatures and pressures and from these we can produce any sort of 
mean or distribution function opacity required. 
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Kurucz, Furcnlid, iirauit, and Testerraan have produced a high- 
resolution, high-signal-to-noise FTS solar flux atlas for the wavelength 
range 296 to 1300 nm. The atlaa is plotted as pseudo-residual flux, and 
a table is given with conversion factors to the Neckel and Labs (Solar 
Physics, 74, 231-249, 1981) absolute calibration. The atlas is being 
published with help from the National Solar Observatory and is being 
sold to recover the cost of printing, A few introductory pages from the 
atlas are attached, but we are able to send only one complete bindery 
sample to NASA headquarters at this time. This atlas has been invaluable 
in testing our line data, 

Larry Petro, Warren Rosen, and Peter Foukal have been searching for 
solar limb-darkening variations that presumably would reflect changes in 
the structure of the solar atmosphere from changes in convective efficiency 
(Petro, Foukal, .Rosen, Kurucz, and Pierce, Ap.J,, 283, 426-438, 1984). 
Kurucz computes the theoretical derivatives for such variations for 
various treatments of convection. Unfortunately our old models do not 
reproduce the solar limb-darkening and do not include the molecular opacity 
that would be the most sensitive component to the variation. Our new model 
does have the molecules but it still does not predict the correct limb- 
darkening. We hope that our next try with the improved atomic opacity 
will finally do the job. 
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RECENT THERMAL MODELS OF THE CHROMOSPHERE 


Eugene Ii, Avrett 

Harvard- Smithsonian Center for Astrophysics 
Cambridge, MA 02138 


ABSTRACT 

This paper describes recent research on the calculation of solar chromos- 
pheric models. The models included here are based on the observed spectrum, 
and on the assumption of hydrostatic equilibrium. The calculations depend on 
realistic solutions of the radiative transfer and statistical equilibrium equations 
for optically thick lines and continue, and on including the effects of large 
numbers of lines throughout the spectrum. Although spectroheliograms show 
that the structure of the chromosphere is highly complex, one- dimensional 
models of particular features are reasonably successful in matching observed 
spectra. There has been considerable recent progress in applying such models 
to the interpretation of chromospheric observations. 

This paper is organized as follows: 

I. The Ca II Resonance Lines 

II. The 40-140 nm Continuum 

III. Model Calculations 

IV. The Temperature Minimum Region 

V. Energy Balance 

VI. Questions and Problems 

The final section lists .a number of unanswered questions and unsolved 
problems that need further study. ' 


I... The. Qa ILRraQDaii.cf-Lliigg 


In the visible part of the solar spectrum a few prominent lines, such as the 
resonance lines of Ca II, have such a high opacity near their central 
wavelengths that the line cores are formed at chromospheric heights where the 
temperature increases in the outward direction. As a consequence, these lines 
show emission near their central wavelengths. Figure 1 shows the observed flux 
profile of the Ca II K and H lines centered at 393.37 nm and 300.85 nm, respec- 
tively. Two emission peaks close to each central wavelength are barely discer- 
nible in these plots of the entire line profiles. The central portion of the K line 
is shown on an enlarged scale in Figure 2, from White and Livingston (1981), 
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Fig. 1. — Ca II K and H flux profiles, obtained with the Kitt Peak Fourier 
Transform Spectrograph (Kurucz el al. 1984). The flux is plotted {vs. 
wavelength in nm) on a linear scale between zero and a pseudo- continuum 
level fitted to the observed flux maxima at 378 and 402 nm. 



INTENSITY 


The lower profile in Figure 2 corresponds to the quiet- Sun intensity at 
disk center, rather than to the flux from the entire disk. The Intensity profile 
from a bright plage region also appears for comparison. 

Solar plages are much brighter than the quiet Sun, but quiet areas have 
brightness variations as well. Figure 3 shows the distribution of intensity pro- 
files of the Ca II II line observed by Cram and Daine (1083). Ten profiles are 
plotted; the lowest corresponds to the faintest 10% and the highest to the 
brightest 10% of the observed quiet region. 



WAVELENGTH 


Fig. 2. — The Ca II I< intensity profile near the line center, from the quiet Sun 
and from plage regions (White and Livingston 1981). 



These results are given in residual intensity units, specifically, tin inien it 
relative to the local continuum level, defined by White and Sue mot o (IlHtK) 
as a smooth curve passing through the estimated continuum values in tin 
far wings of the H and K lines. Paaachoff (1071) studied the relationship 
between these residual units and the measured absolute intensities, and 
concluded that the White and Suemoto continuum intensity is 
(d.5±0.!) X 10 6 ergs cm 2 s'* sr' 1 A '. Tliis value is used to convert the resi- 
dual intensity scale in Figure 3 to absolute units. 

In order to compare with symmetric profiles from static atmosphere calcu- 
lations, the red and violet halves of the line within ±0.7A of the line center 
have been averaged (without including the absorption feature at AX ■■ 0.65 A. 
X wm 306.80 nm, in the violet wing). Figured shows the highest .the median, 
and the lowest profiles from Figure 3, averaged in this way, >n an absoluw 
intensity scale. The corresponding brightness temperatures are shown on tin 
right side of Figure d. 
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Fig 3 — Quiet-Sun H-line profiles observed with high spatial resolution, 
grouped in 10 equal -area bins, from the sequence of profiles ordered according 
to brightness 


The intensity as a function of AX in the line wing beyond the H 2 emission 
peak (at AX ~ O.10A) provides direct information about the temperature as a 
function of optical or geometrical depth, since the intensity successively further 
from the line center is formed at (successively greater depths. The Hj intensity 
minimum is formed at the temperature minimum In the atmosphere; the inten- 
sity rises between Hj and H 2 (at least in the brighter components) because of 
the temperature rise in the low chromosphere. The central reversal occurs 
because the line center is formed high in the chromosphere where the electron 
density is too small to maintain local thermodynamic equilibrium (LTE), and 
photon scattering causes the line- center excitation temperature to be much 
smaller than the local kinetic temperature. 



Fig, 4 . — Smoothed and symmetrized profiles from Figure 3 plotted on an 
absolute intensity scale, with corresponding brightness temperatures shown on 
the the right. 



The brightness temperature at the JI| minimum, or inflection, corresponds 
approximately to the minimum temperature in the atmosphere along the line of 
sight. From Figure 4 this Implies that the average quiet Sun has a minimum 
temperature of about d*100K, and that the faintest 10% and the brightest 10% 
of quiet areas have temperature minima of about 42QUK and -1G00K, respec- 
tively, A more precise interpretation of the H and Iv intensity minima must 
account for the effects of partial frequency redistribution. See Lirtsky (1085) for 
a review of this topic. 

These observations provide information about the temperature minimum 
and low chromosphere because of the high opacity values near the Ca II line 
center. Continuous opacities of the same magnitude occur in the far infrared 
and in the ultraviolet. 

In the far infrared near 150 fim the observed continuum brightness tem- 
perature passes through a minimum value of about ‘MOOIC and becomes larger 
at larger wavelengths. Approximately the same minimum value is observed 
near 100 nm in the ultraviolet, and the observed continuum brightness tem- 
perature becomes larger at smaller wavelengths. The continuum observations 
in the 150 fim and 100 nm regions, as of 1980, have been compiled by Ver- 
nazza, Avrett, and Looser (1070, 1081’, hereafter VAL-I1 and VAL-III), More 
recent, observations are shown <n Figure 22 of the present paper. 

IK The 40- HO nm Continuum 

The lower part of Figure 5 shows the increase in chromospheric brightness 
temperature for extreme ultraviolet wavelengths \<1*10 nm. The upper part of 
this figure shows the corresponding disk- center intensities. Here we see the 
broad wings and central emission of the hydrogen La line, centered at 
121.57 nm. The far wings of the H and IC lines in Figure 1 decrease toward the 
line center because of the outwardly decreasing photospheric temperature, while 
the La- wings increase toward the line center because of the outwardly increas- 
ing chromospheric temperature. The three abrupt intensity increases in Fig- 
ure 5 near 110, 91, and 50 nm are caused by the onset of photoionization from 
the ground levels of carbon, hydrogen, and helium, respectively. 

Apart from the prominent La line, the observations in Figure 5 refer only 
to the continuum values between the many emission lines in the spectrum. The 
complete spectrum, with the lines labeled, appears in Figure 6. These observa- 
tions contain a wealth of information, both about the structure of the solar 
atmosphere and about the detailed radiative and atomic interactions that take 
place. 

Figure 6 shows the spectrum from average quiet regions, from localized 
active regions, and from ’’coronal holes” where the high- temperature coronal 
lines are weak or absent. An off- limb spectrum showing strong coronal lines is 
also included, The count rate at any given wavelength in Figure 6 is propor- 
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tional to intensity. At 00 nm the coronal- hole intensity is about a factor of 2 
less than the quiet-Sun intensity, but active regions are 5 to 6 times brighter 
than quiet regions at this wavelength. 

Different components such as quiet and active regions and coronal holes 
clearly have different intensities. As in the case of the Ca II lines, substantial 
brightness variations also occur within quiet regions. Figure 7 shows the pat- 
tern of brightness variations over 5'x5' quiet areas seen at extreme ultraviolet 
wavelengths. The pattern is a "network” of bright chromospheric emission 
enclosing darker areas often called ’’cells”. 

From histograms of the brightness distribution of such patterns, six com- 
ponents A through F have been selected, ranging from a dark point within a 
cell (component A, the darkest 8%) to a very bright network element (com- 
ponent F, the brightest 4%). At 00.7 nm the darkest 8% of the observed quiet 
area has an intensity about 0.-1 times the average, while the intensity of the 
brightest 4% is almost 3 times the average. See VAL-III for details. 

Figure 8 shows the Intensity distribution for components A through F at 
those wavelengths in the ‘10-1*10 nm range where continuum intensities can be 
easily identified. The component C distribution represents the average quiet 
Sun and is essentially the same as the spatial -average intensity distribution in 
Figure 5. 

III. Model Calculations 


Now consider the determination of the temperature and density variation 
along the line of sight from these and other observations. The observed bright 
ness temperatures at various wavelengths indicate the approximate kinetic tem- 
peratures that occur at unit optical depth for these wavelengths, apart from the 
differences between brightness and kinetic temperatures caused by departures 
from LTE. 

The optical depth at any wavelength depends on the wavelength- 
dependent properties of the absorbing atoms and on the atomic number densi- 
ties along the line of sight. The number densities depend on a balance between 
the solar surface gravity, which compresses the atmosphere, and the gas pres- 
sure, including kinetic and magnetic effects, which expands it. 

In the static case with only gas pressure opposing gravity, the pressure 
varies with geometrical depth x according to the hydrostatic equilibrium solu- 
tion 


p(x) — p(x 0 ) exp 


f dx' 
{H(x') ‘ 


( 1 ) 
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In this case, H is the static pressure scale height H a given by 

H = ( 1+1 'He+ r e)kT 

s_ gm H (l+4r IIe ) 

where r 1Ie is the helium fraction (assumed to be 0.1) and r e — n c /njj is the cal- 
culated electron number density relative to hydrogen. If there are turbulent 
motions, the totai pressure p is the sum of the gas pressure and the turbulent 
pressure 

Pturb = (l/2)pv t 2 urb , (3) 

where p is the gas density and v lurb is an assumed turbulent velocity parameter. 
Then the pressure scale height is 

H = H s +(l/2g)v t u, b . (4) 

Values of v tUrb can be obtained from the nonthermal line broadening velocities 
that must be introduced to account for the observed Doppler widths of various 
lines. These are the only velocities considered here. Magnetic- Held effects are 
not considered. 

Given T, v turb( and r e as functions of x, one can determine p(x) from equa- 
tion (1), apart from a normalizing constant. The electron fraction r e depends 
on T and p and other calculated parameters; the pressure equation is normal- 
ized so that continuum optical depth unity at 500 nm occurs at a fixed photos- 
pheric reference depth, namely at the zero point of the adopted geometrical 
height scale. Thus p(h) can be obtained given T(h) and v turb (h), where h is the 
height above r soo = 1. 

Given p and T, one can compute the number densities, opacities, and 
monochromatic optical depths. The various atomic number densities are 
obtained by excitation and ionization calculations that depend not only on p 
and T but also on the radiation intensity as a function of depth and 
wavelength. The radiation intensity at each depth depends, in turn, on the 
number densities at surrounding depths. These interactions are determined by 
solving the equations of statistical equilibrium and radiative transfer. Thus, 
given T(h) and v tl]rb , one can then calculate all needed atmospheric parameters 
and, finally, the emergent intensity as a function of wavelength. 

The computed chromospheric spectrum is relatively insensitive to the 
chosen values of v turb , even though the geometrical thickness of the atmosphere 
increases substantially as v lurb is increased. Figure 9 shows the adopted v lurb (h) 
distribution (from VAL-III), based on Doppler line widths that are observed to 
exceed thermal widths. For simplicity this distribution is used for all brightness 
components A through F. The atmospheric model then depends only on T(h). 

Figure 10 shows the temperature distribution obtained by trial and error so 
that the computed emergent intensity agrees with various observations of the 



average quiet Sun, The temperature is shown as a function of both h(km) and 
the column mass m(g cm" 2 ), which equals p/g from the hydrostatic equilibrium 
equation. Figure 10 also shows the approximate depth regions where the vari- 
ous continua and line components originate. 

The ultraviolet continuum near 160 nm, the far infrared continuum near 
150 /im, and the minima in the wings of the Ca II and Mg II resonance lines all 
arc formed In the temperature minimum region near h — 500km, The Ca 11 
and Mg II emission peaks, the hydrogen Ha line center, the millimeter contin- 
uum and carbon 110 nm continuum, and the wing of La at 1A and 5A from 
line center all are formed in a broad 5500-7000K temperature plateau that 
characterizes most of the chromosphere. The Ca II and Mg II resonance line 
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Fig. 9. The adopted turbulent velocity distribution compared with values 
given by various authors, 











centers, the centimeter continuum, and the hydrogen Lyman continuum all are 
formed at the base of the chromosphere- corona transition region where the 
temperature increases abruptly. The central part of the La- line and the 3 cm 
continuum are formed at 20,000-30,000K temperatures. 

The plateau near 25,00QK has been introduced to account for l) the Lo 
and L/? integrated intensities, 2) the central absorption observed in these lines, 
and 3) the short- wavelength intensity in the Lyman continuum. The plateau 
could be replaced by a more gradual temperature rise between 10,000 and 
30,000K, adjusted to give the observed values of the integrated Lyman line 
intensities, but the calculated L /? line would not have the central absorption 
that is observed, and La would have less central absorption than observed. As 
discussed later, no physical mechanism has neen identified to explain such a 
plateau. A one- dimensional model may not properly characterize the real 
three-dimensional atmosphere in this respect; the plateau may signify only that 
the line of sight passes through material of the indicated thickness at the indi- 
cated temperatures. 

Figure 10 also shows the depths at which the 11a line core and wing are 
formed. The reason for the separation of the core and wing depths is that the 
hydrogen n — 2 level population is negligible throughout the temperature 
minimum region. Spectroheliograms taken in a narrow wavelength band cen- 
tered on Ha therefore show chromospheric structure and behavior with little 
photospheric contamination. Displacing the band slightly from the line center 
enhances the relative brightness of different chromospheric features having dif- 
ferent velocities along the line of sight. Spectacular photographs of chromos- 
pheric inhomogeneities are obtained in this way. 

The temperature distribution in Figure 10 corresponds to the component C 
observations in Figure 8. The temperature distributions found in the same way 
for components A through F are shown in Figure 11. Note that the six curves 
ail resemble each other in showing a gradual increase in the 6000-7000K tem- 
perature range and a very rapid increase above 8000K. Brighter components 
have higher temperatures at all depths and have chromosphere- corona 
transition regions located closer to the photosphere, at higher pressures. Note, 
however, that the same v turb (h) is used for all six components. 

Active- region models based on the same v turb (h) have temperature 
enhancements relative to Model F that are comparable to the temperature 
differences between Models A and F, and have transition regions located even 
closer to the photosphere. See Basri et al, (1970). However, increasing the 
values of v turb (h) would move the transition region to greater heights. 
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rV J ....T.lifi..To imwrntnr(* Minimum Renlon 

The VAL-I1I temperature- minimum values were based on the 135 -IBS nm 
ultraviolet continuum observations shown in Figure 12. The temperature distri- 
butions for Models A through F in the minimum region are shown in Figure 1J'-. 
Model C is similar to the Harvard -Smithsonian Reference Atmosphere, HSRA, 
of Gingerich el al, (1971), which was based on similar ultraviolet continuum 
observations. Figure 13 also shows the higher temperature- minimum values of 
Ayres and Linsky (1976), determined from a fit to the Ca II and Mg II lines. 
They obtain a quiet-Sun minimum temperature of 4450K compared with the 
Model C value of 4170K. 

As discussed in the VAL-III paper, the 4170K minimum temperature 
appears to be too low for the following reasons. 1) The observed Ca II K] 
minimum intensity is larger than that calculated from Model C by a factor of 
about 1.0. 2) Rast, Kneubuhl, and M tiller (1078) found that the minimum 
observed brightness temperature in the far infrared is 4530(4-100, -150)K near 
130 //m, which is larger than the value computed from Model C by about 200K. 
(This is a higher observed minimum than suggested by the data compiled in 
1976 by VAL-II). 3) Comparisons with the ultraviolet center- to- limb obser- 
vations of Samain (1978, 1080) show that the Mode! C temperature decreases 
too rapidly with height on the photospheric side of the minimum and then 
increases too rapidly on the chromospheric side. Thus the calculated models of 
the temperature minimum region based on the 135-168 nm continuum intensi- 
ties in Figure 12 appear to be incorrect. 

A major uncertainty in the calculations is how to deal properly with the 
many linos throughout the spectrum, The lines are important 1) in calculating 
the photoionization -rate integrals that effect the ionization equilibrium of Si, 
Mg, Fe, C, Al, and other atoms, and 2) in calculating the spectrum for direct 
comparison with observations. 

R. Kurucz has computed new opacity tables based on more than 1.7 X10 7 
atomic and molecular lines. These tables can be used to compute a detailed 
spectrum for comparison with observations, as shown in Figure 14 (from Avrett 
and Kurucz 1083), The darker curve in this figure represents the observed, 
uncalibrated, intensity from the Preliminary Kitt Peak Photoelectric Atlas 
(Brault and Testerman 1972) while the calculated intensity is shown as the 
lighter curve. (See Kurucz and Avrett 1981 for an explanation of the printed 
atomic and molecular identifications.) In the wavelength interval shown in Fig- 
ure 14 the line list appears to give an adequate representation of the observed 
spectrum. Comparisons in other wavelength regions show that many lines in 
the observed spectrum are not included in the current opacity tables, but these 
new opacities give much better results than the approximations used previously. 


I (Wcm ^sr 


ORIGINAL PAG2 S3 

Op POOR QUALITY 



Fig. 12. — Observed and calculated intensities in the 100 nm wavelength ranee. 
From VAL* III, 
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Fig. 14. — Calculated spectrum (lighter curve) in the 380 nm wavelength range 
compared with the spectrum observed by Drault and Testerman (1972). From 
Avrett and Kuruez (1983), 



/Vs reported by Avrott, Kuruez, and Looser (108*1), the VAL-1II calcula- 
tions were repeated using the total line opacity (as a function of temperature* 
and pressure) sampled at 7,000 wavelengths between MO and 005 nm. The 
more realistic line opacities imply higher minimum temperatures than before, 
because l) the ionization equilibrium of SI, Mg, and Fe in the temperature- 
minimum region is closer to LTE than in the earlier calculations, and 2) weak 
lines affect the interpretation of the apparent continuum In the 140-108 nin 
region. 

Figures 15 und 10 show a compilation of observed disk- center intensities, 
and corresponding brightness temperatures, between 1 ftm in the Infrared and 
125 nm in the ultraviolet. Figure 15 Includes two sets of observations, the 
broad -band intensity distribution and the maximum intensities at selected 
wavelengths between absorption lines. 

The crosses in Figures 15 and 10 are the newly- calculated broad- band 
intensities; each cross (70 in all) represents the average of 100 individual inten- 
sities. To compensate for missing lines, all opacities for X <380 nm have been 
multiplied by the factor 1.2, but there are still some discrepancies in the 
210-230 nm region. 

These calculations were made with a new Model C' having a minimum 
temperature distribution dose to that of Ayres and Linsky (1070), as shown in 
Figure 17. This figure also shows the temperature as a function of height for 
new Models A' and FL The three new models have the same 
temperature- height distributions as Models A, C, and F in the chromosphere 
above the minimum region. 

Figure 18 shows a comparison of the II line profiles from Models A', C\ 
and F' with the three observed profiles shown earlier. We attempted to adjust 
each T(h) to fit these observations but, in view of numerous uncertainties, did 
not feel justified in choosing a minimum temperature for Model C' any higher 
than that of Ayres and Linsky, rince they obtained good agreement with the 
observed K line profile. 

These calculated PI j minima generally agree with the corresponding 
observed features to within about 20%, much better than the disagreement by 
the factor 1.0 noted earlier for Model C. Part of the remaining discrepancy 
may be due to uncertainties in the theory of partial frequency redistribution 
applied to the formation of these lines. (See VAL-III and Linsky 1085,) 

The new results are compared with the 160 nm observations in Figure 10, 
which gives the Model C' intensities that should correspond to Samain’s 
observed values between various strong lines (mostly emission lines) in this 
wavelength range. This figure also shows the corresponding intensity distribu- 
tions calculated from Models A' and F'. The new higher- temperature average 
model fits the spatial- average observations as well as before, because the more 
accurate treatment of line opacities forces the ionization equilibrium of silicon 
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Fig. 15. — Observed central Intensities and the corresponding brightness tem- 
peratures between 200 nm and 1 £im, Crosses represent the broad- band calcu- 
lated intensities described in the text, Original figure and references from 
VAL-I1. 
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Fig. 16. — Observed central intensities and the corresponding brightness tem- 
peratures in the 125-200 nm range. Crosses represent the broad-band calcu- 
lated intensities. Original figure and references fom VAL-II. 
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and other atoms to be closer to LTE. In the new calculations, the silicon con- 
tinuum source function at h = 500 km in this ultraviolet wavelength range is 
only about 30% larger than the Planck function, vs. 3 to 4 times larger in the 
VAL-II and -III calculations, and in similar calculations by Cuny (1071), 

The A 1 and F' ultraviolet intensities in this figure are those corresponding 
to the lowest- and highest- decile calcium line profiles in Figure 18, At 168 nm 
the F' intensity is about 3.5 times larger than the A' intensity. Such a differ- 
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Fig, 17. — Temperature distributions in the minimum region for Models A' , 
C, and F' compared with Models A, 0, and F from VAL- III, and with the 
Ayres and Linsky model (shown as a dashed line). 
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ence seems consistent with the quiet -Sun intensity variations with position 
observed by Co^k, Brueckncr, and Bartoe (1983), a.s shown in the upper part of 
Figure 20. The lower part of this figure shows the intensity variations with 
position across a plage region. Similar intensity scans at 160 nm by Foing and 
Bonnet (198*1) are shown in Figure 21 . Much could be learned from a study of 
the correlations between such high -resolution intensity variations in the ultra- 
violet minimum region and simultaneous high- resolution intensity variations in 
the Ca II K t or Mj minima. 



Fig, 18. — H-line profiles calculated from Models A 1 , C', and F' (lighter 
curves) compared with the Cram and Dame observed profiles (darker curves). 
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Fig, 19, — Ultraviolet intensities calculated from Models A 1 , C 1 , and F* (with 
lines included throughout the spectrum as described In the text) compared with 
observations and earlier calculations. 



Finally, the calculated and observed brightness temperatures in the far 
infrared are compared in Figure 22. This figure shows the observed minimum 
value 4530(+100, -150)K at 130 /<m from Rast, Kneubuhl, and Muller (1078) 
and recent observations at 50, 80, and 200 ^im reported by Degiacomi (108-1). 
(See also Degiacomi ei al. 108-1.) The open circles represent quiet- Sun observa- 
tions, while the filled circles represent the higher values obtained when active 
regions were present within the area observed. Individual (small-scale) active 
regions cannot be resolved with the 00 cm telescope used for these submillime- 
ter wavelength observations. 

Model C' appears to be in adequate agreement with the quiet- Sun obser- 
vations in Figure 22. Thus, a minimum temperature of 4400 or 4500K for the 
average quiet Sun seems consistent with the Ca II, 160 nm, and 150 (im data. 
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Fig, 20. — The absolute intensity and corresponding brightness temperature at 
168,6 nm from quiet and plage regions at two times, 30 s apart, from slit spec- 
tra obtained during a rocket flight of the Naval Research Laboratory’s High 
Resolution Telescope and Spectrograph (Cook, Brueckncr, and Bartoe 1983). 
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Fig. 21. — Intensity scans at 100 nm across an average quiet-Sun area and an 
active region. The intensity scale is in arbitrary units. The 2-sigrna level due 
to photographic noise is indicated on the left. From Foing and Bonnet (1084). 
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Fig. 22. — Far-infared brightness temperatures calculated from Models A', C' , 
and F' compared with the earlier model calculations, with the observed 
minimum value 4530 (-f 100, -150) K at 130 (im from Rast, Kneubiihl, Muller 
(1978), and with the observations of Dcgiacomi (1984) at 50, 80, and 200 /im. 
(Open circles represent quiet- Sun observations; filled circles represent observa- 
tions with active regions included.) 
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The two preceding sections have summarized how hydrostatic-equilibrium 
models of the chromosphere and temperature minimum region are determined 
from various observations. Such models should provide Information about chro- 
mospheric heating. 

Given the depth variation of both the number densities of the atomic and 
molecular energy states, and the radiation intensities in the important lines and 
continua, one can calculate the net radiative cooling rates 4> u ; and <I> m in ergs 
cm -3 s" 1 as functions of depth for lines between levels u and l , and for con- 
tinua associated with levels m. The sum of all of these cooling rates must be 
the non- radiative energy input at each depth required to maintain the given 
temperature distribution, assuming a steady state. In radiative equilibrium 
without any non- radiative energy input, the sum o.f these cooling rates would 
have to be zero, but the chromosphere is clearly not in radiative equilibrium, 

The net radiative cooling rate, or radiative flux divergence, for continua is 
given by 

= MM Siz-Jyjdt' , (5) 

m 

where k„ is the continuous opacity (in cm" 1 ), S„ is the source function, and is 
the mean intensity of radiation. For a line, equation (5) becomes 

* u / — \\u [n u (A u j + B u , J) - n,B /u jJ, (0) 

where n y and n / are the upper and lower level number densities, where A u / , 
B u / , and B/ u are the Einstein coefficients for spontaneous emission, stimulated 
emission, and absorption, and where J is the mean intensity weighted by the 
line opacity and integrated over the line according to 

J = (f K^di^/if K u du ) . (7) 

Using the Model A', C', and F' atmospheric parameters, the net radiative 
cooling rates have been calculated for l) the lines and continua of atomic 
hydrogen, 2) the bound- free and free-free continua of H“, 3) the resonance and 
subordinate lines of Ca II, 4) those lines of Mg II, and 5) the lines and continua 
of He, He II, Ca, Mg, Si, Si II, C, C II, Al, Fe, Na, K, and O. These five groups 
are called H, PI", Ca II, Mg II, and “other contributions,” respectively. 

Figure 23, consisting of the plots shown on the following 18 pages, shows, 

first for Model C 1 and then for Models A 1 and F', the following six plots: the 


Fig. 23. (next 18 pages) — Net radiative cooling rates, first for Model C' and 
then for Models A' and F', showing the total <E> with the five contributions to 
the total, followed by five plots showing the details of these contributions. 
(The symbol BC on the hydrogen plot refers to the Balmer continuum; the 
integer pairs in the Ca II and Mg II plots refer to upper and lower line indices.) 
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total <I> with the five contributions to the total, and then details of each of the 
five contributing groups. The upper and lower panels on each page show posi- 
tive and negative values of <I>, on logarithmic scales, plotted against height in 
km above continuum optical depth unity at 500 nm. 

The total rate shows three Important characteristics: a) the maximum near 
800 km in the low chromosphere and the gradual decrease of 4>< fola | above this 
height, b) the maximum In the lower part of the chromosphere -corona transi- 
tion region at temperatures between 20,000 and 3O,O00K, and c) the negative 
values of ^Total near ‘100 km l* 1 the temperature minimum region. (The large 
positive values below the temperature minimum region have less significance; 
they occur because <I> Is proportional to the rapidly increasing number density 
in this region and because the photospheric temperature distribution has not 
been adjusted to obtain ^Total ~ In Model C' the total hydrogen number 
density at 100 km is 2500 times larger than the value at 1000 km.) 

The maximum of <I> Tola | in the low chromosphere and its gradual decrease 
with height shows how mechanical energy is dissipated with height to maintain 
the chromospheric temperature rise. Theories of wave energy dissipation must 
be consistent with such calculated values, as discussed in the papers by Stein 
and others In this volume. The chromospheric mechanical energy input escapes 
primarily as Ca II and Mg II line emission. 

The maximum of <h Tola | In tlm transition region corresponds primarily to 
emission in the hydrogen La line. The area under this maximum is 
~3X10 5 ergs cm -2 s“ l in the ease of Model C 1 , which agrees with the observed 
integrated La intensity in quiet regions near the disk center. This La feature 
seems required whether there Is a temperature plateau in the transition region 
or a more gradual temperature rise instead (or a succession of interfaces 
between chromospheric and coronal temperatures intersected by the line of 
sight), Rabin and Moore (1084) have recently proposed a localized heating 
mechanism to account for this required energy. Alternatively, one might sup- 
pose is that the energy 3X10 5 ergs cm" 2 s" 1 is dissipated in the corona and is 
then conducted downward to the lower transition region to be radiated away by 
La. A difficulty with such a mechanism is that the required conductive energy 
transport can take place only if the temperature gradient is very large (assum- 
ing classical conduction theory), while a much smaller gradient is needed in 
order to produce the observed emission. This difficulty may be resolved if 
non-Maxwellian effects are properly included, as discussed by Shoub (1982, 
1983) and Owocki, Canfield, and McClymont (1985). 

The negative values of $ Tola) in the temperature minimum region are per- 
plexing because negatives imply non-radiativc energy extraction, i.e,, refrigera- 
tion. Possible explanations include 1) the effect of time-dependent wave 
motions that, on the average, indeed remove energy from the temperature 
minimum region, 2) the possible omission of positive contributions that would 


make the total rate positive, and S) the need to introduce inhomogeneities that 
have a more extreme range of variation than components A 1 through F'. These 
three possibilities are discussed below. 

1) If waves travel through the atmosphere and the temperature at each 
position varies as s function of time, the time* averaged computed spectrum 
may have characteristics of .Vigil and low temperatures that cannot Le matched 
by an average temperature, or by separate components such as A' through FV 
Recent calculations by Muchmore and (Jhnschneider (1984) indicate that the 
negative values of ^xotai in the temperature minimum region might be 
explained in this way. 

2) There is a positive <h at 400 km from the group, of 13 “other contribu- 
tors” Me through 0, but the negative <I> from at this height has an absolute 
value typically 10 times larger than the sum of all positive values. However, 
potentially important contributors not included in this group are a) prominent 
lines of molecules such as carbon monoxide, and b) large numbers of weak 
atomic and molecular lines. 

To estimate these effects, Kurucz’s line- opacity tables sampled at 1055 
representative wavelengths between 144 nm and 10 // m have been used to com- 
pute the monochromatic flux divergence 

<t>v = 47r« l XS„-J„) (8) 

as a function of depth. The integral of <p u over all u should give <f» Totat . This <f> u 
is integrated from 144 nm to 10 /.tm, thus excluding contributions from 
X<144 nm (e.g., Let) that affect 4>'p otia ] in the middle and upper chromospheie. 
Also, the central features of the Ca II and Mg II resonance lines are not among 
the sampled wavelengths, so that $ Tola] computed in this way does not include 
the terms that are most important in the chromosphere. However, this pro- 
cedure should give a valid result in the temperature minimum region. Fig- 
ure 24 (for Model C') shows ^xotai obtained by Integrating <£„ from 144 nm to 
10 //m. This figure also shows the five component integrals over the 
wavelengths 144-300 nm, 300 nm-l/nn, 1-4.2/um, 4.2-5. 4 /tm, and 
5.4 - 10 /tm. Note that 'h'poial ' s negative at 400 km, just as in Figure 23. 
These negative values are caused mainly by the integral from 300 nm to 1 /t m, 
which is dominated by the H" bound -free continuum. Positive values are 
found for all of the \>1 /tm contributions, including the H" free- free contin- 
uum. These positive and negative values can be understood from the tempera- 
ture dependence of the Planck function in the two wavelength ranges. Fig- 
ure 25 shows (for Model C') the Planck function B, the continuum source func- 
tion S, and the corresponding mean intensity J at 5 /tm and at 500 nm plotted 
against height and monochromatic optica! depth. At both wavelengths, J at 
400 km is about 0.3 times S at h = 0, but S > J at 5 /tm while S < J at 
500 nm, because of the greater change of B with temperature at 500 nm. Thus 
<p u from equation (8) is positive at 5 /tm and negative at 500 nm. The 5 /tm 
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Fig. 25. — Planck function, source function, mean intensity, and optical depth 
at continuum wavelengths 5 fim and 500 nm. 








graph shows continuum data, but the corresponding results for a strong carbon 
monoxide line near 5 /<m are similar, except that k in equation (8) is larger, unit 
optical depth occurs higher, and J is closer to S. The computed emergent 
intensity values in the 4.2 - 5.4 band are shown in Figure 20, in order to 
indicate that the sampled wavelengths include strong carbon monoxide lines. 
These calculations are intended to show that even when the CO lines and many 
other lines are included, the net radiative cooling rate for Model C' remains 
negative in the temperature minimum region. Calculations with Models A' and 
F' lead to the same result. 
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Fig. 26. — Computed intensity values in the 4.2 - 5.4 /x m band, indicating the 
lines of 



3) It may be possible to avoid the negative values of 4> Tola | if the atmo- 
sphere consists of inhomogeneous regions having a greater temperature range 
than components A 1 through F'. Consider a bifurcated model consisting of a) 
background regions having minimum temperatures well below those of 
Model A' , and b) bright points with minimum temperatures much higher than 
in Model F'. For temperatures much less than 4000K the CO lines can become 
strong enough to balance the negative cooling rate due to H^, and a 
radiative- equilibrium solution can be obtained with ^Toia) — 0* For tempera- 
tures near 5000K the negatives due to are reduced in the minimum region 
as S approaches J, and a second radiative- equilibrium solution can be obtained. 
These solutions are discussed by Ayres (1081), Kneer (1983), and Muchmore 
and Ulmschneider (1085), Negative values of <I> ToU | would be avoided if the 
atmosphere consists of these two components, each having 4> Tota | ” 0* Such an 
explanation is supported by the observations of Ayres and Testerman (1981) 
who measured the fundamental vibration- rotation band of CO near 4.6 ftm. 
They found strong CO lines with central brightness temperatures as low as 
3700K near the limb. The substantial limb darkening of these lines was noted 
earlier by Noyes and I-lall (1972). Recent observations by Deming ct al, (1985) 
also provide support for a bifurcated model. They find horizontal thermal fluc- 
tuations of order ±800K in the upper photosphere. 

While a bifurcated model may explain the CO and OH line observations 
and the values of 4> Tola ], such a model seems inconsistent with the Ca II, 
far- infrared, and ultraviolet observations. The H line profiles in Figures 3 and 
4 indicate that the faintest decile of the observed quiet Sun has an average 
minimum temperature of about 4200K, and that the brightest decile has an 
average minimum temperature of about 4 BOOK. The 160 nm observations in 
Figure 20 show similar results. In the far infrared near 150 /um, the observed 
quiet- Sun minimum brightness of — 4400K means that the sum of the bright- 
ness temperatures of each fractional component weighted by its fractional area 
has to be — 4400K, since LTE prevails and since the intensity is proportional to 
temperature in this wavelength region. The 4400K average could be the result 
of half the area at a temperature of 5000K and half at 3800K, for example, or 
the result of 70% at 5000K and 30% at 3000K. On the other hand, 
high- resolution Call and 160 nm observations suggest that most of the 
observed quiet areas have minimum temperatures within the intermediate range 
4200-4G00K. Thus it seems difficult to reconcile a simple combination of 
high- temperature bright points and a low- temperature background with the 
4400K spatial average, the observed brightness distribution of the Ca II 
minimum features, and the observed brightness distribution at 160 nm. 

Of the three possible explanations of the negative values of 4> Tota | discussed 
above, the first one also involves hot and cool fractions, which change with time 
and with vertical as well as horizontal position. As discussed by Muchmore and 



Ulmschneider (1085), wave rarefaction can produce localized cooler regions 
which can be cooled further by CO radiation, while compression can produce 
localized hot regions. Their time- dependent solutions of the mass-, momen- 
tum-, and energy conservation equations confirm that the temperature 
minimum region has a bi- stable character, as discussed by Ayres (1981) and 
Kneer (1983). 

While Ca II and 100 nm observations with sufficient spatial resolution may 
rule out a simple bifurcated model with temperature extremes sufficient to give 
«I> Tota l = 0 in both components, such observations may be consistent with a 
dynamical interpretation, for which the temperature varies with time and with 
vertical position, and for which hot and cool elements can occur alternately 
along the line of sight. 

VI. Questions and Problems 


I conclude this paper with a list of unanswered questions, unsolved prob- 
lems, and unfinished work, related to the particular topics discussed here. 

1) Detailed chromospheric models, based on observations (and the assump- 
tion of hydrostatic equilibrium), have been given for various components of the 
quiet Sun. Similar detailed calculations are needed to improve the available 
models of active regions (Basri et al, 1979, Orrall 1981), coronal holes, sunspots 
(Avrett 1981, Lites and Skumanich 1982, Maltby et al. 1984), and flares 
(Machado et al, 1980, Ricchiazzi and Canfield 1983, Canfield, Gunkler, and 
Ricchiazzi 1984). 

2) Cook, Brueckner, and Bartoe (1983) have demonstrated that the chro- 
mosphere is extended to greater heights than hydrostatic equilibrium implies, 
and that the active -region chromosphere extends much higher than the quiet 
chromosphere. What causes the scale height to be greater in active regions? 
How do magnetic fields affect the structure of quiet and active regions? (See 
Spruit 1981; Stenflo 1983.) 

3) Only small portions of the available ultraviolet spectral data (such as in 
Figure 6) have been used to determine detailed atmospheric models. Further 
systematic work is needed to study the observed optically thick lines that are 
formed in the chromosphere and lower transition region (Lites, Shine, and Chip- 
man 1978; Tripp, Athay, and Peterson 1978; Lemaire and Gouttebroze 1983; 
Roussel- Dupre 1983b). Further work also is needed to improve the theory of 
line formation with partial frequency redistribution. See Roussel-Dupre (1983a) 
and Linsky (1985). 

4) The source of energy that results in La emission has not yet been iden- 
tified. Is this energy first deposited in the corona, and then brought to the base 
of the transition region by non- Maxwellian effects (Shoub 1982, 1983; Owocki, 
Canfield, and McClymont 1985), or does some mechanism deposit the energy 
directly in the La region (Rabin and Moore 1985)? 


5) Ultraviolet temperature- minimum observations made with high spatial 
resolution, as in Figures 20 and 21, show a brightness distribution similar to 
that obtained from high- resolution observations of the Ca II minimum features. 
Simultaneous observations of these spectra would provide useful constraints on 
the spatial distribution of hot and cool components of quiet and active regions. 

6) In addition to disk observations in the far infrared, observations at the 
limb, and particularly eclipse oceultations at wavelengths > 150 /mi, provide 
valuable information about the structure of the temperature minimum region 
and low chromosphere. See Lindsey et al. (1983, 1084), and Cram (198-1), 
These far- infrared observations also show that the chromosphere is more 
extended than the hydrostatic equilibrium solutions indicate, 

7) Regarding the negative values of the net radiative cooling rate <J> To tal * m 
the temperature minimum region for Models A' through F', it has been sug- 
gested that higher temperatures would give a radiative equilibrium solution 
with $ Total = 0 (due to reduced H£} heating), -and that lower temperatures 
W'ould give a second ^Total = 0 solution (due to enhanced CO cooling). These 
two radiative equilibrium solutions need to be calculated accurately, including 
the effects of line blanketing. See Kalkofen (1985). 

8) The Ca II, 160 nm, and 150 (x m observations may rule out a simple 
bifurcated model composed of high -temperature bright points and a 
low -temperature background, with $ Total = 0 for both components. On the 
other hand, these observations may not rule out a dynamical model having 
similar properties. These ideas need further study. 

I am grateful to Robert Kurucz and Rudolf Loeser for their contributions 
reported here, and to Bruce Lites and Robert Noyes for their comments on the 
manuscript. This research was supported by NASA Grant NSG-7054. 
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SOLAR FLUX ATLAS FROM 296 TO 1300 NM 


X. INTRODUCTION 

Thin atlas proaonto the aolar flux spootrum observed with high roso- 
lutlon aij't high algnal-to-noiso using tho Fourier transform spootromotor 
(FTS) on the HoMath Solar Toleaoopo at Kltt Peak National Observatory. 

Tho FT'S allowed ua to greatly lmprovo tho nignal-to-noiso and wavelength 
coverage compared to tho earlier high roaolutlon atlas produced by Bookers, 
Bridgos, and Gilliam (1970 using a grating apootrograph at Sacramento Peak. 
Information about the' observations is summarized in Table 1, Eight • ->onn, 
odd-numbared 1 to 15, wore mado of the integrated solar disk with ovw it .pping 
bandpasses that oover tho spootrum from 295 to 1300 nm. High points from 
these scans are plottod in Figures l and 2 for froquonoy and wavelength. 

At tho same time scans, oven-numborod 2 to 16 , wore made for a standard lamp. 
Tho ultraviolet soana have rosolving powor of approximately 348000 and tho 
infrared scans approximately 522000. Peak signal-to-noise for each scan 
ranges from 2600 to 9000 as listed in Table 1, At a particular wavelength 
tho noise can be best estimated from tho spootrum itself. 


II. WAVELENGTH REDUCTION 

The frequency scale for each 3can must be normalized by a multipli- 
cative factor near unity that depends on the setup and alignment of 
the FTS. Ideally this factor would be determined by observing through 
a gas absorption oell for a gas with well-known lines in the bandpass. 

Wo determined this faotor from the spectrum itself by using the 0„ line 
at 688.38335 nm and defining its air wavelength to be that given By Pierce 
and Breokenridge (1973), The spectra were converted from frequency to air 
wavelengths using Edleu's (1966) expression. As thi3 line appears only in 
scans 11 and 13, it could not bo used for the other scans. Scan 15 was 
set by aligning terrestrial lines common to soan 13. Since the other 
soana did not have suitable terrestrial lines, they were corrected in the 
solar rest frame by aligning relatively clean solar linos in overlapping 
scans. The solar rest frame was dftfcerminod from radial velocities com- 
puted with a program that we obtained from Mark Reid of the Center for 
Astrophysios. It has errors on the order of 4 m/a beoause it does not 
inolude planetary perturbations. The radial velocities for the beginning, 
end, and middle of each 3can are listed in Table 1. The c.iange in radial 
velocity while a scan was being made resulted in about 0.2 km/a broadening 
of the solar lines. Wo used the middle value in shifting the spectra. 
Operationally, the frequency correction factor wa3 al3o treated as a 
Doppler shift and specified in km/3. These values are listed In Table 1. 
Considering all the fitting and shifting, our final wavelengths in scan 1 
may have errors as large as 0.1 km/s. If the wavelength of any one line 
in each 3can can be determined more accurately, the whole 3can can be 
rescaled to that accuracy. 


III. REDUCTION TO RESIDUAL FLUX 

To produce a residual flux spectrum we made a zero-point correction. 
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divided out the standard loop, fitted a pseudo-continuum to oaoh soan, 
and than wo combined the scans together into a continuous whole. 

Hon-linoiritiea in the detector introduce oystematio errors in the 
strong. low-i v <aquenoy componenta of the Fourier transform that produoe a 
varying zero-point error in the spectrum. Thia error can be avoided by 
using a narrow bandpass to reduce the signal. We corrected the spectrum 
by fitting lie ,ar segments to opaque terrestrial line centers that were 
wido enough that ringing was not a problem. Ringing is dearly visible 
in the expanded scale plots in Figure 7 for terrestrial lines that are 
not resolved. Only scans 11, 13, and 15 had terrestrial lines that 
allowed this correction, and for those scans the maximum corrections were 
0.45%, 0.25%, and 0.90%, respectively. Ho estimate the residual error 
to be about 0.10% . Scans 1 to 9 were not corrected and we cannot rule 
out the possibility of 0.25% zero-point errors. In particular, the cores 
of the Ca II H and K lines are uncertain by this amount. 

In conjunction with each 30lar scan, a standard lamp was scanned 
Hith the spectrometer in the same configuration. We divided the solar 
scans by the lamp scans to remove the effects of any small features in 
the bandpass limiting filters that were usee with each scan. 

Figure 3 illustrates our method for subjectively fitting a pseudo- 
continuum to eaoh scan. For each scan the highest point out of every 100 
points wa3 plotted and a curving pseudo-continuum was subjectively fitted 
to the highest high points. Each scan was normalized to its pseudo-con- 
tinuum so that it became residual flux with the continuum level at 1.0. 
Note in Figure 3 that the pseudo-continuum is obviously far from being 
the real continuum foP scans 1 to 7 . 

Figure 4 illustrates the ( subjective fitting together of the 3cans 
into a continuous whole. We made plots of the overlapping scans and chose 
a high region where the scans coincided as the crossover point. The wave- 
lengths where the cuts were made are 329.897, 378.2914, 401.965, 473.8, 
576.5, 753.9, 999.7 nm. These plots also give an indication of repro- 
ducibility and noise level but these are regions near half-maxima of the 
bandpasses where the bandpasses have steep slopes of opposite sign. The 
quality of the data at the centers of the bandpasses is much higher. 
Ideally measurements like these should be taken with rectangular band- 
passes that overlap in order to avoid these problems. 

There are systematic differences between scans. First, zero-point 
corrections could be made only on scans that had completely opaque ter- 
restrial lines that are supposed to produce zero signal. No zero-point 
correction could be made for the blue and ultraviolet scans and we cannot 
rule out ex’rors as large as 0.25% . If neighboring scans have different, 
non-negligible errors the residual fluxes are inconsistent across the 
join. Second, these data are of such high quality that very weak terres- 
trial features can be significant. Neighboring 3can.s were taker* through 
different atmospheric conditions with different amounts of water vapor 
and ozone. For example, at 402.07, 402.14, and 999.9 nm the spectra are 
actually different for this reason. Finally, the merged spectrum has 
solar wavelengths so terrestrial lines are Doppler shifted from eaoh 
other in neighboring scans by the amounts listed In Table I. 
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The final residual flux spectrum is plotted twiuo, in Figure 5 at 
40 nm per page, and in Figure 6 at 4 nm per page in the visible and at 
8 nm per page in the infrared. In Figure 6 the top 10% is also plotted 
at an expanded scale that shows weak features. 


IV. IRRADIANCE CALIBRATION 

The calibration is based on the irradianoo derived by Neckel and 
Labs (1984). Their data are listed in Table 7 . and shown in Figure 7. 
Their bandpasses provide continuous coverage from 329 to 666 nm and 
sporadic points out to 1247 nm. From the atlas we computed the mean 
value of residual flux R for each bandpass assuming that the bandpasses 
were exaotly rectangular. The width and R for each bandpass are listed 
in Table 2. The actual bandpasses were slightly rounded and have 3ome 
wavelength error, but errors from these effaots are Inconsequential. 

We divided the observed irradianbe by the mean R to find the irradl- 
ance pseudo-oontlnuum level listed in Table 2. 

Our value of R inoludes the terrestrial lines. Neckel and Labs' 
corrections for terrestrial line absorption in their observations are 
also listed in Table 2. As they chose their bandpasses in the red and 
Infrared to have few or no terrestrial absorption lines, the correc- 
tions vary from 0 to 2% and are typically on the order of 0.1% . We 
plan to determine this for our scans by synthesizing both the solar 
and the terrestrial spectrum, but we have not yet done so. We have 
assumod that our terrestrial absorption is approximately the same as 
that of Neckel and Labs' and that ours Is negligible when theirs is 
negligible. Figure 7 includes a plot of the pseudo-continuum values 
corrected for thi3 absorption. Ideally this should be a smooth curve, 
but, in fact, it has structure at the several per cent level. 

The dip at 600 nm is probably caused by our poor treatment of the 
ozone absorption there, which, by accident, falls between scans 9 and 
11. The ozone feature is the Chqppuis bands which run from 450 to 750 
nm. The cross-section is tabulated by Vigroux (1953) and is roughly 
plotted in a^figure in Goody (1964). The maximum absorption coefficient 
Is 5.E-21 cm which should produce about 5% absorption at 600 nm in 
our atlas. Our signal-to-noi3e is high enough that we should be able 
to see the structure. Vigroux also gives data between 407 and 447 nm. 
Apparently there are no measurements of ozone cross-sections between 350 
and 400 nm and longword of 750 nm. At our signal-to-noise, ozone may 
be significant at all wavelengths. We would very much like to see a 
thorough, high-resolution laboratory study. 

In the infrared we subjectively fitted an expression for the 
pseudo-continuum to the measured points for the range 746.5 to 1300 nm, 

IRRADIANCE=0.80933E18/W**5/ 

(EXP(1.4384E7/W/(6892.094-2.691159*W+.001669325*W**2))-1.) 
for wavelength W in nm and irradlance in pW/cm 2 /nra. 

Below 330 nm we pieced together a calibration using Broadfoot (1972) 
whose long wavelength limit is 319.9 nm and Arvesen et al. (1969) who 


start at 300 nm and go up. Wo plotted Broadfoot, Arvesen, and Neckel 
including overlaps and oonoluded that the Arvesen wavelength scale was 
shifted by +0.4 nm and that their irradiance should be reduced by about 
10 per Qont to be in agreement with Neckel and Lab3 a-ii > 'th Broadfoot, 

We made Neokel-like 2.1 nm bands for Broadfoot and 2 nm t>. : ds for Arvesen 
et al. These are listed in Table 2 and plotted in Figure 7. We U3ed 
Broadfoot for 296 to 310 nm and Arvesen ot al. multiplied by 0.9 from 
312 to 32t nm. We do not trust this ultraviolet calibration to bettor 
than 25 per oent, both from the observations and from the very poor 
pseudo-continuum fit in this region as 3hown in Figure 3. 

To oonvert the atlas from pseudo-residual flux to irradiance, we 
computed a pseudo-continuum value, either from the fit in the infrared 
longward of 746.5 nm, or from linearly interpolated bandpass values, 
and then multiplied by the pseudo-residual flux. Table 3 lists the 
calibration of the pseudo-oontinuura level of the atlas in |iW/cm Z /nm 
for every 0.05 nm in the visible and every 0.10 nm in the infrared. 

One |iW is 10 ergs/s. For flux in erg3/om z /s/nm at the surface of the 
sun multiply the irradiance by 10.+(1.495985E13/6.9598E10)**2 - 462020. 
For astrophysical flux F divide the flux by n. .For flux moment H divide 
the flux by 4n. 


V. MAGNETIC TAPE LISTING 

A magnetio tape listing of the atlas may be obtained by sending a 
2400 foot tape to Dr. Robert Kurucz, Center for Astrophysics, 60 Garden 
Street, Cambridge, Massachusetts 02138. The tape will be written 9 track, 
1600 bpl blocked ASCII with 80 characters per record and 50 records per 
block. The first file will contain wavelength, residual flux pairs. The 
second file will be a program that can read the tape and oonvert it to 
irradiance if desired. The third file will be the residual flux inter- 
polated to uniform point spacing in wavelength. 
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Table 

1. Record of 

the Observations 



scan 1 

scan 3 

3can 5 

3can 7 

scan 9 

scan 11 

scan 13 


29 6. - 329. 

897 - 37 8„ 

2914 - 401 

.965 - 473 

8 - 576.5 - 753. 

9 

i 

99 

29113.63238 

25147.01194 

24649.74115 

17171.55918 

12485.51810 

8624.80162 

7103.72630 

i 

.03645036 

.03528395 

.02979210 

.05516111 

.02381424 

.01815228 

.01421564 

i 

i 

695994 

696954 

696362 

680994 

1042714 

1045794 

1045974 

3 

f 

348000 

348480 

348184 

340500 

521360 

522900 

522990 

1 

5000 

9000 

3000 

3000 

3000 

3000 

2600 

4 

1 

06/22/81 

06/21/81 

06/22/81 

11/23/80 

03/24/81 

03/25/81 

03/25/81 

; 

f: 

V 

42205.78 

43752.23 

34183.89 

38243 .02 

45035.48 

38297.73 

47186.79 


49637.39 

54066.75 

41652.48 

44168.07 

52805.08 

45284.97 

54248.59 

i 

.01039 

.01439 

.01039 

.0079 8 

.00719 

.00639 

.00639 

1 

26 

36 

26 

20 

18 

16 

16 

1 

UV 

UV 

UV 

silver 

silver 

silver 

silver 

■s 

Si diodes 

Si diodes 

Si diodes 

Si diodes 

Si diodes 

Si diodes 

Si diodes 

r-t 

Cor2985+ 

CuS04+ 

CS7-51+CS4-97 

CS7-59+CS4-97+ 

GG435+ 

VP700(15)+ 

RGN9 

f 

CuS04(10) 

CS7-54 

WG345 

GG375 

CS4-96 

OG570 


§ 

$ 

-2365.44 

-806.28 

-10386.13 

-5406.42 

208.00 

-6512.91 

2050.36 

i 

4737.39 

9179.02 

-2918,65 

189 .79 

7651.57 

148.14 

9113.65 


.409003 

.409102 

.409018 

-.324822 

.028530 

.034821 ' 

.035529 

| 

.228211 

.165580 

.700693 

.998589 

.529622 

.702180 

.541084 

l 

.338631 

.606078 

.256191 

.915586 

.744536 

.522780 

.8129 82 

1 

1.02 

1.01 

1.30 

1.84 

1.15 

1.30 

1.16 

1 

1.05 

1.21 

1.03 

1.63 

1.35 

1.15 

1.45 

S 

601.3 

601.5 

601.3 

599 .7 

601.3 

602.2 

602.2 

1 

i 

0.029 

0.076 

-0.154 

-0.476 

0.492 

0.309 

0.557 

1 

0.222 

0.330 

0.015 

-0.320 

0.704 

0.502 

0.741 

1 

0.125 

0.203 

-0.070 

-0.398 

0.598 

0.406 

0 . 649 

l 

-0.75 

-0.75 

-0.75 

0.25 

0.15 

0.30 

0.30 
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of the Observations 


scan 11 

scan 13 

scan 15 


753 

.9 - 999 

.7 - 1300. 

adopted wavelength limits (nm air) 

8624.80162 

7103.72630 

5247.91886 

wavenumber of first point 

.01815228 

.01421564 

.01154261 

dispersion (cm-1) 

1045794 

1045974 

1044994 

number of points in 3can 

522900 

522990 

522500 

spectral resolving power 

3000 

2600 

3000 

max signal-to-noise 

03/25/81 

03/25/81 

05/11/81 

date 

38297.73 

47186.79 

37338.11 

starting time (seconds) 

45284.97 

54248.59 

45506.13 

stopping time (seconds) 

.00639 

.00639 

.00719 

tlme/polnt (seconds) 

16 

16 

18 

number of scans 

silver 

silver 

silver 

beamsplitter 

Si diodes 

Si diodes 

InSb 

detector 

VP700(15)+ 

RGN9 

GaAs+ 

optical filters 

OG570 


CS5-57 


-6512.91 

2050.36 

-6895.55 

hour angle at start (seconds) 

148.14 

9113.65 

944.94 

hour angle at end (seconds) 

.034821 ' 

.035529 

.314112 

mean declination (rad) 

.702180 

.541084 

.532358 

zenith distance at start (rad) 

.522780 

.812982 

.251236 

zenith distance at end (rad) 

1.30 

1.16 

1.16 

airmass at start 

1.15 

1.45 

1.03 

airmass at end 

602.2 

602.2 

599.4 

atmospheric pressure (torr) 

0.309 

0.557 

0.184 

radial velocity at start (km/s) 

0.502 

0.741 

0.397 

radial velocity at end (km/s) 

0.406 

0.649 

0.291 

radial velocity mean (km/s) 

0.30 

0.30 

0.00 

frequency correction (km/s) 


I 
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